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Abstract 

Using the [SU(3)i J x S , C r (3)R,] g iobai X [SU (3) y] local chiral Lagrangian with hidden local symmetry, 
we evaluate the cross sections for the absorption of eta meson (rj) by pion (71"), rho (p), omega 
(u), kaon (K), and kaon star (K*) in the tree- level approximation. With empirical masses and 
coupling constants as well as reasonable values for the cutoff parameter in the form factors at 
interaction vertices, we find that most cross sections are less than 1 mb, except the reactions 
prj — * KK*(KK*), lot] — * KK* (KK*), K*n — > pK, and K*n — * ujK, which are a few mb, and the 
reactions ttw — > KK and Kn — ► irK, which are more than 10 mb. Including these reactions in a 
kinetic model based on a schematic hydrodynamic description of relativistic heavy ion collisions, 
we find that the abundance of eta mesons likely reaches chemical equilibrium with other hadrons 
in nuclear collisions at the Relativistic Heavy Ion Collider. 

PACS numbers: 25.75Nq, 12.39.Fe, 13.75.Lb,14,40.Aq 
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I. INTRODUCTION 



Besides the intrinsic interest of studying eta meson production in heavy ion collisions 
[lj], knowledge of the eta meson dynamics in the hot dense matter formed in these colli- 
sions is also important for understanding other observables measured in experiments. For 
example, to extract information on the rho meson in-medium properties from the dilepton 
spectrum it emits requires reliable knowledge on the number of eta mesons produced in 
heavy ion collisions as eta mesons through their Dalitz decays contribute significantly to 
low-mass dileptons measured in these collisions Q]. Also, the two-pion correlation function 
measured via the Hanbury-Brown-Twiss interferometry is affected by pions from eta decays 
P). Understanding the eta meson dynamics is thus essential for extracting the size of the 
pion emission source from the pion interferometry. In heavy ion collisions at energies avail- 
able from the Heavy Ion Syncrotron (SIS) at the German Heavy Ion Research Center (GSI), 
whose dynamics is dominated by baryons, the final eta number is mainly determined by its 
interaction with the nucleon through the N(1535) resonance j^, which has a branching ratio 
of about 35% decaying into pion and nucleon. At higher energies from the Super Proton 
Syncrotron (SPS) at the European Organization for Nuclear Research (CERN), the matter 
becomes more dominated by mesons, and the eta-meson interaction is thus important in 
determining the final eta yield |5J. For heavy ion collisions at the Relativistic Heavy Ion 
Collider (RHIC) at Brookhaven National Laboratory (BNL), the final hadronic matter is 
largely made of mesons, and it is even more crucial to have a good knowledge on the cross 
sections for eta meson absorption by mesons. 

Since there is no empirical information on these cross sections, theoretical models are 
needed to evaluate their values. A possible model is the [SU (3)l x 5 , f^(3)R] g i bai x [SU (3) v] local 
chiral Lagrangian with hidden local symmetry [6j. This model has been used to study 



hadronic matter In this paper, we shall use this model to evaluate the cross sections 
for the absorption of r\ meson by 7r meson to the final states of KK, KK*(KK*), K*K*, 
pcu, and irp; by p meson to the final states of KK, KK*(KK*), K*K*, pp, ttu, and tttt; 
by ijj meson to the final states of KK, KK*(KK*), K*K*, and up; and by K and K* 
to the final states of K(K*)tt, K(K*)p, and K(K*)u. In evaluating the cross sections for 
these reactions, we do not include the effects due to scalar mesons a(980) and a(1450) as 
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studies based on the non-linear chiral Lagrangian have shown that contributions from these 
resonances are unimportant compared to those from the vector meson exchanges and contact 

n 

interactions |9(. These cross sections are then used in a kinetic model to study their effects 
on eta meson abundance in heavy ion collisions at RHIC. 

This paper is organized as follows. In Section |H1 we give a brief description of the 
[S77(3)l x S"?7(3)R]giobai x [SU ( 3) v] local chiral Lagrangian with hidden local symmetry. The 
interaction Lagrangians that are relevant for describing the absorption of 77 meson by n, 
p, u, K, and K* are derived in Section IIHI The cross sections for these reactions are 
then evaluated in Section IIV1 In Section the time evolution of rj meson abundance in 
relativistic heavy ion collisions is studied in a schematic model. Finally, a summary is given 
in Section IV II Explicit expressions for the amplitudes for 77 meson absorption are given in 
Appendices lAllEl 



II. CHIRAL LAGRANGIAN WITH HIDDEN LOCAL SYMMETRY 

Based on the hidden local gauge symmetry, Bando [6.J has constructed an effective theory 
for hadrons interacting at low energies. In this approach, vector mesons are introduced as the 
gauge bosons of the hidden local symmetry of the nonlinear chiral Lagrangian. For hadrons in 
the SU(3) multiplets, which are relevant for present study, it is constructed with two SU(3)- 
matrix valued variables and £,r(x) that transform as £,l,r(x) — * r( x ) = h(x)^L,R g\ r 
under h(x) G [SU(3)v]i oca j and g^R G [SU(3)L,R] g iobai- The resulting chirally invariant 
Lagrangian is 

£ = £ A + a£ v -^(F, u Fn, (1) 

with 

£a = -\f:((D^ L -d-D^R-Q 2 

= d^Vy - d v Vp - ig[V^ V v \. (2) 

In the above, f n is the pion decay constant, the symbol (■■) denotes the trace of 3 x 3 matrix, 
and = dp — igV^ is the covariant derivative, with the dynamical gauge bosons of the 
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hidden local symmetry identified with the nonet of vector mesons, i.e., 



/ 1 n° 4- 1 i i 



V 



p + K* + \ 



if* 



(3) 



The above effective Lagrangian with a = 2 is known to give the universality of vector meson 
coupling, the Kawarabayashi-Suzuki-Riazuddin-Fayyazuddin (KSRF) relations JlOf . and the 
vector meson dominance of the pseudoscalar meson electromagnetic form factor 
Fixing the [SU(3) v }i oca \ gauge by 



i\ = £ R = £ = exp(z$// 7r ] 



(4) 



where $ is the nonet of pseudoscalar Goldstone bosons 

if 



\/2 



V 



7T 



71 n 



if 

■■jsi + y/lJ] 



(5) 



leads to the usual lowest-order chiral Lagrangian 

/(f) c . t\ _ ; y 



^ ((^ • e f - ^u f ■ o 2 ) = f -f (dj/pu*) , 



(6) 



-20 c 



HQ 



with [/ = £ 2 . In obtaining Eq.©, we have used the empirical mixing angle 9 
between the octet r) 8 and the singlet % to obtain the physical r\ and rj via 77 ps (2y / 2^8 + ^i)/3 
and V « + 2 v / 2?? 1 )/3. 

The SU(3) symmetry breaking effects are taken into account by introducing in the La- 
grangian the mass term 



with the mass matrix M. given by 

M = diag(m 2 , m 2 , 2m 2 K — m 2 ). 
Also, La and Cy are modified as follows jl^j ]: 

£ A + ac a = - l -fi ((d^ l . - ^ • d) 2 a + + e^el) 
l v + A£ y = -1/2 ud^ l ■ ei + ■ ^) 2 (i + ex^ei + e^ei) 



(7) 



(8) 



(9) 



with 6a(v) — diag(0, 0, ca(v)), where ca(v) are re& l parameters to be determined by empirical 
hadron masses and decay constants. 

Expanding the pseudoscalar fields in the first equation of Eq.® shows that the kinetic 
terms for K and rj are renormalized by the symmetry breaking term. The proper kinetic 
terms can be recovered by rescaling the K and 77 fields according to [14f 



VI + c A K -> K and y 1 + -c A rj -> 77. (10) 
As a result, the kaon and eta decay constants are related to the pion decay constant by 



f K = VTTca'U and f v = yl + ^cu/w, (11) 
and masses of vector mesons become different, i.e., 

2 2 m lf* m | /-2 2 / 10 \ 

m„ = m = = — — aizq . 12 

p w l + c v l + 2c v J * y K J 

Expanding the Lagrangian up to four meson fields leads to interaction Lagrangians of 

the VPP, VVV, WPP, and PPPP types, i.e., 

Cvpp = -^({[<M M $],^}(l + 2e y )>, 

Cvvpp = ~ag 2 (V^iev® 2 + 2$e y $ + &e v ) 
Cpppp = ^(d^d^$-d^d^<£ 2 ) + ^-(M¥ 

-- i (^$^$(2^$ + $ 2 e A + e A $ 2 )) . (13) 

J n 

In the chiral Lagrangian with hidden local symmetry, vector mesons have also been 
included through the anomalous interaction of the VVP type in order to take into account 
the breaking of local chiral symmetry Q|. Including further the flavor breaking through 
a term £lCwz£b + £rCwz£l with e wz = diag(0, 0, c wz ) [14], the total anomalous interaction 
Lagrangian then reads as 

C VVP + ACwp = 2g VVP e^(d,V„(l + 2e wz )d x V a $). (14) 



III. INTERACTION LAGRANGIANS 



Inserting Eqs.tjUJ) and © in Eqs.(JT3J) and (fUj) and rescaling the K and 77 meson fields 
according to Eq. fjlOj) . we obtain following interaction Lagrangian densities that are relevant 
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7] meson absorption by n, p, u, K, and K*\ 

ag 



£pKK 
£wKK 

£k*Kti 

r 



£pK*K v 
£-u>K*Kr\ 

£>-KT]KK 



-PPV 



■'UIUIT) 



£pK*K 
£wK*K 
£>K*K*tt 
£>K*K*ri 



2 

ag 1 



= %- 



p» ■ (vr x 9 M ?f), 

(Krd^K - d^KrK) ■ 



4 1 + CA 

ag 1 



= % 



4 1 + CA 

] + Cv K*t ■ (Kd^-d^K^+R.c, 
4 V 1 + c a 

i% • 1 + * V 2 M ( Kd ^ - d^Krj) + H.c, 



+K*»(T-f?d»K;-T-d,frK;)}, 

C^k*k* = t^[(d,K^K;-K^d,K;)^ + (K^d,uj u -d,K^uj u )K^ 



ag 



cv 



2V / 6 V / (l + c A )(l + |c A ) 



ag 



cv 



2V6 X /(1 + CA )(1 + | CA ) 



— . 1 + ){Krd tx K -8^ + d^KrK ■ d^nrj) 

+Kt8 ix K ■ Trd^rj + d^KrK ■ nd^r] - (2 + 3c A ) KrK ■ d^rj 
- 2d li Kf& t K ■ ttt] + m^KrK ■ nrj 



gppr] 



^d,p v ■ d X p a T), 



gp^ a(3Xa d^p u dxuj a ■ 7T, 

gg£g 7 =L — (l^ • r9 A K; + d,p u ■ d x KlfK) , 
y J: e^8^ v (Kd x K* a + d x K* a K) , 



2^ *K*7] C wz 

V6 y'l + |c A 
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In the above, r are Pauli matrices for isospin; and tt and p denote the pion and rho me- 
son isospin triplet, respectively; and K = (K + ,K°) T and K* = (K* + ,K*°) T denote the 
pseudoscalar and vector strange meson isospin doublet, respectively. 

For the parameters in the interaction Lagrangians given in Eq. (|15Jl . we choose a = 2 to 
recover the vector dominance model and the KSFR relation. Using the empirical value f n = 
92.4 MeV for pion decay constant, m p = 776 MeV for rho meson mass, and nix* = 892MeV 
for K* mass, the KSFR relation (Eq. 1)12)1 ) then gives a a universal vector coupling constant 
of g = 5.94 and the parameter cy = 0.32. We note that a slightly larger value of cy = 0.36 
was obtained in Ref.jSj by fitting the empirical phi meson mass = 1020 MeV. Since only 
the K* is needed in present study, we choose to have the correct K* instead phi meson mass. 
From Eq . 1)11)1 . the parameter ca = 0.49 is obtained from the empirical kaon decay constant 
fx = 1-22/,. 

It is interesting to know that above parameters lead to a width T p = g 2 {m 2 p — 
4m^) 3//2 /(487rm^) ~ 147 MeV for rho meson decay to two pions and a width Yk* = 
[g\l + c v ) 2 )/{l + c A )]{[m 2 K , - {m K + m n ) 2 ][m 2 K „ - {m K - m n ) 2 ]} 3 / 2 /(16nm 5 K .) ~ 49 MeV 
for K* decay to Kit, which are close to the empirical values of T p = 150 MeV and Tk* = 51 
MeV. 

For the anomalous coupling constants Qyyp, they are related to the universal vector 
coupling g and pseudoscalar decay constant by [lQj . 

9vvp = q ?- , P = n, r], K (16) 

87T Jp 

and from the ratio between experimental decay widths of K*° K°j and K*^ K^j, 
the parameter c wz = —0.1 has been determined jl^ . 

IV. ETA ABSORPTION BY MESONS 

Since the temperature of the hot hadronic matter produced at RHIC is above T = 125 
MeV, it consists of not only the lightest it mesons but also heavier p and u mesons as well 
as the strange mesons K and K*. Eta mesons can thus be absorbed by all these mesons. 
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A. Born diagrams 



Diagrams for these reactions are shown in Fig. ^for rj absorption by it to the the final 
states of KK, KK*(KK*), K*K*, pu, and 7tp; in Fig. |2]for rj absorption by p meson to 
the final states KK, KK*(KK*), K*K*, pp, -nuo, and 7T7t; in Fig. El for 77 absorption by uo 
meson to the final states KK, KK* , K*K* , and up; and in Figs. 0] and [S] for rj absorption 
by K and K* to the final states of irK, irK*, pK, pK*, uK, and uK*. 



K K K K K K 




I I 

TC r\ TZ T\ 

C4a) (4b) 



7T P 
P 

7r r\ 

C5a) 

FIG. 1: Diagrams for rj absorption by tt meson. 

Using the interaction Lagrangian densities given in Section IIII1 we have derived the 
amplitudes for all tree-level diagrams shown in Figs. Q|SJ In general, the amplitude for a 
process n is given by 

M n = A\£M%- x ^e kXk ...ea v (17) 

where % runs through all the subprocesses in each reaction, and denotes the polarization 
vector of external vector meson j. The factor A is either a matrix element rf- of the Pauli 
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FIG. 2: Diagrams for r] absorption by p meson. 

matrices or the kronecker delta 5 a b or the antisymmetric tensor ie a bc- It takes into account 
the isospin states of the particles in a reaction, with a, b, and c denoting those of isospin 
triplet Ti and p meson, and i and j those of isospin doublet K and K* . Explicit expressions 
for these amplitudes are given in Appendix [X] for rj absorption by ix meson, Appendix El 
for 7] absorption by p meson, Appendix O for 77 absorption by to meson, Appendix O for rj 
absorption by K meson, and Appendix |E] for 77 absorption by K* meson. 



B. form factors 



To obtain the full amplitudes for these reactions, one needs in principle to carry out a 
coupled-channel calculation in order to avoid the violation of unitarity. Such an approach 
is, however, beyond the scope of present work. To prevent the artificial growth of the tree- 
level amplitudes with the energy, we introduce instead form factors at interaction vertices 
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FIG. 3: Diagrams for 77 absorption by u meson. 

and treat their cutoff parameters as parameters. Specifically, the form factors are taken to 
have the same forms as used previously in studying J/ip absorption 16, ^| and charmed 
meson scattering by mesons, charmed meson production from photon- and proton-proton 
reactions Q], pentaquark baryon production from photon- and hadron-proton reactions 
|2fll |. and strangeness-exchange reactions between mesons and baryons j^jj. For three-point 
vertices, i.e, p-mr, pKK, uKK, K*Kn, K*Krj, ppp, pK*K*, uK*K*, pun, pprj, uur], pK*K, 
u)K*K, K*K*7r, and K*K*rj in the t— and w-channel processes, the form factors are taken 
to have the form 

Fs(q) = a4^' (18) 

where q 2 , taken in the center of mass, is the squared three momentum transfer. For four- 
point contact vertices, i.e., pK*Kr], uK*Kr], and irr]KK, they are taken to have the form 



A? \( A3 



1 



(19) 



A? + (q 2 >; VAi + (q 2 )/' 
where Ax and A 2 are two different cutoff parameters at the three-point vertices present in 
the processes with the same initial and final particles, and (q 2 ) is the average value of the 
two squared three momenta at the three-point vertices. 

For three-point vertices in s-channel processes, we introduce, however, a covariant form 
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FIG. 4: Diagrams for 77 absorption by K meson. 

factor, i.e., 

where s is the square of the center-of-mass energy and m is the mass of the intermediate- 
state particle. This form factor ensures that it would not affect the predicted decay widths 
of p and K*, which are close to the empirical ones, as s = m 2 in these processes. We note 
that unrealistic large cross sections for s— channel processes at high energies can also be 
prevented by taking into account the vacuum widths for all channels that are open at a 
gl ven . a s in hJJ on phi m e S on interactions in hot hadronic n,atte, 

As in previous studies on hadronic reactions, we use A = 1 GeV in calculating the cross 
sections for eta absorption by mesons. Since the value for A might depend on the reactions, 
we also study how our results are affected by changes in its value. 
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FIG. 5: Diagrams for r/ absorption by K* meson. 
C. cross sections for eta absorption by mesons 

The isospin- and spin-averaged differential cross sections for above reactions are given by 

da„ 1 



dt 64nsp?NiNi 



-l-Mnl 2 , (21) 



where |.M n | 2 denotes the squared amplitude obtained from summing over the isospins and 
spins of both initial and final particles, and can be evaluated using the software package 
FORM 22]. The factors iV 7 = {2I X + 1)(2J 2 + 1) and N s = {2S 1 + 1)(2S 2 + 1) in the 
denominator are due to averaging over the isospins l\ and I 2 as well as the spins Si and 5*2 
of initial particles, while is their 3-momentum in the center-of-mass frame. 

Integrating the four momentum transfer t leads to the following total cross sections: 

a n = -J—Pl\M n \ 2 , (22) 
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where pf is the 3-momentum of final particles in their center-of-mass frame and \M n \ 2 is 



related to \M„\ 2 in Ea.(l2Tl) bv 



\M n \ 2 = ^ / dn\M n \*F\ (23) 



with F denoting the appropriate form factors at interaction vertices. 




FIG. 6: Cross sections for eta absorption by ir (left panel), p (middle panel), and u (right panel) 
mesons as functions of total center-of-mass energy s 1//2 above the threshold s^ 2 of a reaction. 

In Figs. Uland [?[ we show, respectively, the cross sections for 77 absorption by tt, p, and uj 

mesons and by K and K* as functions of total center-of-mass energy s 1//2 above the threshold 

1/2 

s of a reaction. Aside near the threshold of a reaction, where the cross section can be very 
large or small depending on whether the reaction is exothermic or endothermic, most cross 
sections are less than 1 mb, except the reactions prj — > KK*(KK*), oorj — > KK*(KK*), 
K*rj — > pK, and K*r] — > uK, which are a few mb, and the reactions irrj — > KK and 
Krj —>■ ttK, which are more than 10 mb. The large cross sections for the reactions tti] —>■ KK 
and Kr] — > nK are due to the presence of the nr]KK contact interaction in their amplitudes. 
Values of the cross sections depend, however, on the value of the cutoff parameter in the 
form factors. We find that increasing the value of A to 2 GeV increases the cross sections 
by about a factor of two, but decreasing its value to 0.5 GeV reduces the cross sections by 
about a factor of four. 
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TABLE I: Values of parameters a, b, c, d, and e in Ea.(|24[) for parameterizing the matrix elements 
denned in Ea.CT. 



Reactions 


(7 


b 


n 


d 




TTT] - 


+ KK 


16.047 


-468.921 


-3.867 


598 598 


-3.215 


TTT] KK*(KK*) 


5.813 


-4.232 


-0.136 


-1.580 


-5.357 


TTT] — > 


K*K* 


0.115 


7.261 


-13.200 


8.468 


-4.530 


-/rr/ - 


-> puj 


1.539 


2.577 


-1.513 


0.555 


-15.636 


7T77 ■ 


-* Tip 


0.721 


-7.012 


-3.948 


6.293 


-4.323 


pr/ - 


* KK 


2 706 


— 7 532 


—5 657 


4 827 

1 -±. \j£t I 


—4 857 


prj -» KK*(KK*) 


36.139 


20.087 


-18.884 


-17.458 


-2.338 


pq -> 


K*K* 


0.181 


-4.174 


-19.494 


13.793 


-4.941 


PV - 


-> pp 


14.265 


—2.296 


—2.812 


—10.394 


—2.812 


PV - 


-> 7T7T 


0.001 


—0.529 


—8.081 


0.529 


—7.476 


PV - 


-> 7TW 


2.007 


-0.699 


-5.410 


-0.657 


0.228 


LOT] - 


-> 


1.863 


-0.948 


-6.925 


-0.878 


0.056 


cur] -> KK*(KK*) 


12.179 


7.285 


-19.514 


-5.766 


-2.153 


LOT] — > 


K*K* 


0.060 


-1.434 


-19.684 


4.669 


-4.948 


LOT] ■ 


— > 7T/J 


7.784 


-0.675 


-5.356 


-6.369 


0.0317 


Kr] - 


-» vrK 


— 33 196 


78 96^ 


—0 460 


— 1 5 97^ 


— 1 373 


Kr] - 


-> 7rK* 


00425 


-18.683 


-5.715 


18.688 


-5.511 


Krj - 


-► pK 


6.819 


-5.863 


-0.871 


-0.955 


1.118 


Kr] - 


-» pK* 


0.288 


-0.413 


-15.160 


1.593 


-5.047 


Kr] - 


-» WK 


3.124 


-2.091 


-0.936 


-1.032 


0.556 


Kr] - 


-» 


-0.389 


-0.322 


-1.785 


1.182 


-0.625 


K*r] 


-> VrK 


2.078 







—2.077 


-2.998 


K*r] - 


-» vrK* 


0.078 


9.298 


-4.638 


0.514 


-36.521 


K*r] 


-> pK 


34.849 


10.952 


-15.036 


-16.308 


-2.838 


K*r] - 


-> pK* 


2.999 


-1.694 


-27.299 


1.181 


-5.701 


K*r] 


-» 0)K 


11.591 


3.936 


-15.140 


-5.276 


-2.807 


K*r] - 




0.966 


-0.522 


-27.134 


0.409 


-5.647 
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FIG. 7: Cross sections for eta absorption by K (left panel) and K* (right panel) mesons as functions 
of total center-of-mass energy s 1 / 2 above the threshold s^ 2 of a reaction. 

The results shown in Figs. E] and [7| can be conveniently reproduced by parameterizing 
the squared matrix elements shown in Eq. (j23j) according to 

\M n \ 2 = a + b(^^j +d (^) mb ' (24) 

The dimensionless parameters a, b, c, d, and e for the different reactions considered in 
present study are given in Table |U 



V. ETA PRODUCTION IN RELATIVISTIC HEAVY ION COLLISIONS 

In this section, we study the effect of eta absorption reactions studied in previous sec- 
tions on the time evolution of its abundance in heavy ion collisions at RHIC. We describe 
the collision dynamics in a schematic hydrodynamic model and assume that the hadronic 
matter is in thermal equilibrium throughout the collision. Furthermore, we assume that 
pions, rho mesons, omega mesons, and kaons as well as their resonances are in chemical 
equilibrium during the collision because of their larger interaction cross sections. For eta 
mesons, we consider two scenarios for their abundance in the beginning of the hadronic 
stage after hadronization of the quark-gluon plasma, i.e., they are either taken to be in 
chemical equilibrium with other hadrons or completely absent. Their abundance during the 
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evolution of the subsequent hadronic matter is then described by a kinetic model that takes 
into account their absorption and regeneration. We also compare the final eta number in 
the two scenarios with that from assuming that eta mesons remain in chemical equilibrium 
throughout the evolution of the hadronic matter. 

A. rate equation 

The density n„ of eta mesons changes in time according to the following rate equation 

□ 

<9>X) = ( 25 ) 

where = 7(1, v) is the four velocity of the hadronic matter fluid element of velocity v 
and 7 denotes the Lorentz factor. The source term \1/ is given by 

* = - X] ( a ari^bcv)n v n a + (<Tbc-+a V v)n b n c , (26) 

a,b,c a,b,c 

where n a , rib, and n c are densities of meson types a, b, and c, respectively. The thermal av- 
eraged eta absorption and production cross sections are denoted by (o~ av ^bcv) and (cr6 C ^ a7? t>), 
respectively, with v the relative velocity of initial two interacting particles. The thermal 
averaged eta production cross sections are related to those of eta absorption cross sections 
by 

(a a ^ hc v)n^ = (<7 bc ^ av v)nt q n?, (27) 

where n eq denotes the equilibrium density, i.e., 

n^= d -f^K 2 {m/T). (28) 

In the above, d and m are the degeneracy and mass of a hadron, T is the temperature of 
the hadronic matter, and K2 is the modified Bessel function of the second kind. Using the 
above relation, the rate equation can be written as 

M (nX) = -£ (a ari ^ bc v)n^(n v - rff). (29) 

a,b,c 
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B. thermal averaged eta absorption cross sections 



With particle momenta in the hadronic matter approximated by the Boltzmann distri- 
butions, the thermal averaged cross sections can be expressed as 23j 

(av) = [4alK 2 (a 1 )a 2 2 K 2 (a 2 )}- 1 

dz[z 2 - (ai + a 2 ) 2 ][z 2 - (ai - a 2 ) 2 ]K l (z)a(s = z 2 T 2 ), (30) 

with ai = rrii/T, z = max(ai + a 2 , a% + 0:4), and K\ being the modified Bessel function of 
the first kind. 



--0 




-pri-KK_ —- pti-KK or h 

pr| — K K pn— ► Jnt 

- pr| — Jtto pTi->pp 



- 1.5 

A 

> 

V 1.0 



0.5 



0.0 l 



n 1 r 



-cup KK 
••«[-» KK or T 
-isq-> itp 



_i 1 1 1 i_ 



10 0.15 0.20 0.25 0.30 0.10 0.15 0.20 0.25 0.30 0.10 0.15 0.20 0.25 0.30 



T(GeV) 



T(GeV) 



T(GeV) 



FIG. 8: Thermal averaged rj absorption cross sections by pion (left panel), rho (middle panel) and 
omega (right panel) mesons as functions of temperature. 



The resulting temperature dependence of (av) is shown in Fig. |H1 for the reactions of 
in], prj, and ujt] and in Fig. El for the reactions of Kr] and K*rj. It is seen that the thermal 
averaged cross sections for most reactions increase with increasing temperature. As for 
cross sections, the reactions Tirj — > KK and Krj — > ttK have largest thermal averaged cross 
sections. 



C. collision dynamics at RHIC 



Since the particle distribution in central heavy ion collisions at RHIC is approximately 
uniform in midrapidity and the geometry of the collision is cylindrically symmetric, it is 
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FIG. 9: Thermal averaged r\ absorption cross sections by K (left panel) and K* (right panel) as 
functions of temperature. 

convenient to use the cylindrical coordinates r, ip, r, and r\ with the latter given by 

T = y/W=#, v =L n l±l. (31) 
2 t — z 

Assuming longitudinal boost invariance and allowing for radial transverse expansion, then 
one has u v — — 0. For uniform density distribution in the transverse plane, averaging 
over the radial coordinate gives 

—^^(Ti^r)^)) = - £ Ka-^XK - O (32) 
rn [T) or a b ^ c 

In the above, R{r) is the transverse radius of the system and (u T ) is the averaged r compo- 
nent of the four velocity and is given by 

( M -) = —_/ drru T {r). (33) 



R 2 {r) Jo 

At midrapidity, the four velocity of hadronic fluid element -u M can be expressed in terms 
of the radial flow velocity f3 r as 

" T= *=7i=* (34) 

With the usual ansatz for the radial velocity, i.e., 

Afr.O = (35) 
18 



we have 



(u T ) = j 1 dy , 1 . (36) 

T - {dR/drfy 



o 



To determine the time evolution of the transverse radius of the fireball, we follow the 
model used in Ref. 



24j . Since we are interested in the time evolution of the eta abundance 
during the hadronic phase of the collision, we start from the end of the mixed phase th and 
write 

R(t) = R h + v h (t - Th ) + ^(t- t h ) 2 . (37) 

In the above, Rh ~ 9 fm and vh ~ 0.4c are, respectively, the transverse radius and flow 
velocity of the fireball at th = 7.5 fm/c, while a = 0.02c 2 /fm is the acceleration in the 
transverse expansion. Values of these parameters are determined from fitting the measured 
transverse energy ~ 788 GeV as well as the extracted freeze out temperature Tp = 125 MeV 
and transverse flow velocity ~ 0.65c of midrapidity hadrons in central Au+Au collisions 
at v /sjvtv = 200 GeV. Assuming that the hadronic matter expands isentropically, the time 



dependence of the temperature of the fireball obtained in Ref. 



24j can be parameterized as 



T(r) =T C - (T H - Tp) f^^L) °' 8 , (38) 

\Tp - T H J 

where Th is the temperature of the hadronic matter at the end of the mixed phase and 
is thus the same as the critical temperature Tq for the quark-gluon plasma to hadronic 
matter transition. As in Ref.j^j], we take Th — Tc = 175 MeV. The freeze out temperature 
Tp = 125 MeV then leads to a freeze out time Tp w 17.3 fm/c. 

D. time evolution of the eta abundance 

To study how the eta meson abundance evolves in time, we consider two scenarios for 
the initial eta meson number, i.e., no eta meson is present at th or the eta meson is in 
chemical equilibrium with other hadrons at th as in the statistical model 2^|. In the first 
case, subsequent hadronic interactions increase the final eta meson number to about 30 at 
freeze out, as shown by the dotted line in the left panel of Fig. El This scenario is, however, 
unrealistic as we expect eta mesons to be appreciably produced during hadronization of the 
quark-gluon plasma. For example, in the quark coalescence model 26], the number of eta 
mesons produced at hadronization ranges from 12 for a small eta meson root mean square 
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FIG. 10: Time dependence of the abundance of midrapidity rj mesons in the hot hadronic gas 
formed from central Au+Au collisions at ^/s~nn = 200 GeV at RHIC for cutoff parameter A = 
1 GeV (left panel), 2 GeV (middle panel), and 0.5 GeV (right panel). Solid and dotted lines 
correspond, respectively, to eta mesons that are chemically equilibrated or absent at the beginning 
of the hadronic phase, while dashed lines correspond to eta mesons that are always in chemical 
equilibrium. 

radius of 0.5 fm to 37 for a larger radius of 1 fm. In the second scenario of equilibrated 
eta mesons, the initial eta meson number is N(th) ~ 35. Including subsequent hadronic 
absorption and regeneration does not change much the eta abundance, and its number at 
freeze out is decreased slightly to about 34 as shown by the solid line in the left panel of 
Fig. This result is consistent with the assumption of the statistical model that hadron 
abundances in heavy ion collisions at RHIC are largely determined at hadronization. For 
comparison, we have also shown in the left panel of Fig. El by the dashed line results 
from the assumption that eta mesons are always in chemical equilibrium with other hadrons 
during evolution of the hadronic matter. In this case, the eta number shows a small initial 
increase that is followed by a decrease to about 32 at freeze out. The similarity among the 
final eta abundance in all three scenarios indicates that the yield of eta meson in relativistic 
heavy ion collisions is not very sensitive to its production mechanism. This is not surprising 
as the cross sections for the reactions ni] — > KK and Krj — > ttK obtained with a cutoff 
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parameter A = 1 GeV in the form factor are sufficient large, implying that eta mesons are 
strongly interacting and thus likely reach both thermal and chemical equilibrium in the hot 
hadronic matter produced at relativistic heavy ion collisions. For a larger cutoff parameter 
of A = 2 GeV in the form factor, these cross sections will be even larger, and we would 
expect an even more similar final eta abundance for the three scenarios, and this is indeed 
seen in the middle panel of Fig. El Even for a smaller cutoff parameter of A = 0.5 GeV, 
which leads to a much smaller cross sections for eta absorption and production, final state 
interactions in the hadronic matter can still keep eta mesons close to chemical equilibrium 
if they are initially in chemical equilibrium as shown in the right panel of Fig. E3 On the 
other hand, the final abundance of eta mesons will be significantly below the equilibrium 
value if they are initially absent in the hadronic matter, which is, however, unlikely as we 
have commented in the above. 



VI. SUMMARY 



Knowledge of eta meson interactions in hadronic matter is not only of interest in its 
own right but also important for extracting information on the properties of the hot dense 
matter formed in high energy heavy ion collisions. Since there is no empirical information 
on the absorption cross sections of the eta meson with pion, rho and omega mesons as well 
as kaon and its resonance, which are the most abundant particles in high energy heavy ion 
collisions, we have evaluated these cross sections based on the tree-level diagrams from the 
[S77(3)l x SU (3 )r] g i bai x [SU(3)y]i oca \ chiral Lagrangian with hidden local symmetry and 
including symmetry breaking effects. Using empirical hadron masses and coupling constants 
as well as reasonable values for the cutoff parameters in the form factors at interaction 
vertices, we find that although cross sections for most eta absorption reactions by mesons 
are less than 1 mb, the reactions tti] —* KK and Krj — > ttK are more than 10 mb. To see 
the effects of these reactions on the yield of eta mesons in relativistic heavy ion collisions, 
we have solved a kinetic equation, based on a schematic model for the dynamics of heavy 
ion collisions, to follow the time evolution of the abundance of eta mesons. We find that 
the final abundance of eta mesons at freeze out of the hadronic matter is close to chemical 
equilibrium, irrespective of the initial value at hadronization of the quark-gluon plasma. 
In particular, if eta mesons are initially in chemical equilibrium with other hadrons, their 
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number is not strongly affected by their subsequent interactions during the expansion of the 
hadronic matter. 
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APPENDIX A: ETA ABSORPTION BY PION 

Using subscripts 1 and 2 to denote the initial-state particles and 3 and 4 for the final-state 
particles in the order from left to right in all the Feynman diagrams shown in Figs. WS^ 

for the polarization vector of vector mesons, and also the usual Mandelstam variables 
s — (Pi +P2) 2 , t = (pi — P3) 2 and u = (pi — P4) 2 , the amplitudes for the absorption by mesons 
can be written explicitly as given below. For propagators, we do not include the width of 
exchanged particles as its effect is negligible due to the large threshold of the reactions. We 
note these expressions are obtained with a = 2 in the chiral Lagrangian. 

1. TTT] — > KK 

The amplitude for this reaction is given by 



M 



KK = TfJM la + M u + M lc ), 



(Al) 



with 



M 



a 




x -g^ u + 



(pi -p?,y{pi -p 3 ) 



(P2 + P4) 




(1 + Cy) 2 




1 



x -g" u + 



(gl -P^iPl -P4) 

m 2 K « 



V 



(P2 + Ps) 



22 



u * "Jit {i f c A jy i -t- ^c A 
+P2 ■ P3 + P2 ■ Pa + (2 + 3c A ) Pi-P2 + 2P3 -pt + ml 



2. ttt) KK*(KK*) 



The amplitude for this reaction is given by 



M 



irri—fKK* (KK*) ~ 1 ij 



T°(MZ a + M% b )e^ 



with 



M 1 



2a 



2b 



99K*K*r, (1 + C V )C V 



X 



Vq ^(i + CA )(i + | CA ) 

(Pi-Ps)Api -Pa)v 



(Pi+PsY 



m 



K* 



t — m 2 K » 



e av '^p,p{p2-pi\ 



9 9k* k*tt 

(1 + Cy)(l +2c wz ) 
" V6 y/(l+C A )(l + lc A ) 



xe a %(pi-p 4 ) a 



-9w> + 



(pi -Pa)v{Pi-Pa)i 

777 * 



u — m K » 



-(P2+P3) 1 



3. TT7] -► if* if* 



The amplitude for this reaction is given by 



with 



MT 



3a 



36 



AC 



(jl {l + Cv? -(2 Pl - P3 r'(2 P2 - P4 y 



Ve(i + c A )yJi + lc A 

9 2 (1 + cv) 2 
V6(l + c A )y/l + lc A 

9K*K*n9K*K* V (1 + 2C WZ )C WZ ^ 



(2pi-p 4 r- 



t — 

1 



u — m 



■(2 P 2-p 3 r, 



K 



X 



V6 



— g 77 ' + 



1 + fcA 



P3a(i?3 -Pl)/3 



(Pi -Ps) 7 (Pi -Ps) 



7' 



._~ 2 € aWp,I/ p4p>(p2 - P4)a', 



9k*k*-k9k*k*t) (1 + 2c wz )c wz o . 
^ ^/l + |c 
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X 



—g~n' + 



(gi -P4) 7 (Pi zl^L 

2 



u — m 



-6 a '^' P3a ,(p2-p3)p>. (A6) 



A'* 



4. 7r?7 — > 



The amplitude for this reaction is given by 



(AT) 



with 



l 4a 



X 



v/l + §<u 

(Pi - P3) 7 (Pl - Ps) 



Q ^ 7 P3a(P3-Pl)/3 



-<?77' + 



t — m? p 



e aV/3 '>4/3'(P2-p4)a', 



r L 1= ^p 4/ ,(p 4 -p 1 )a 



x 



-<?77' + 



(Pi - P4) 7 (Pl - P4)" 



u — mi 



e aV/3 ' 7 W(P2-P3)/J'. (A8) 



5. 7T?7 — > 7T/9 



The amplitude for this reaction is given by 



M 



(A9) 



with 



M£ = 



99ppij 



(pi+Ps)" 



X 



-SW' + 



(Pi -PsMPi -PsV 



i — m 2 , 



e a ^P4 /3 (P2-p 4 ) a . (A10) 



APPENDIX B: ETA ABSORPTION BY RHO MESON 



1. pr] A\K" 



The amplitude for this reaction is given by 



M^ KR = t«(MI + M» h + M^)e lfl , 



(Bl) 
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Mr, 
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M 11 
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2. pr/ — > KK*(KK*) 



The amplitude for this reaction is given by 

■M p rj ^KK*{KK*) = t*(AC + AC + AC + AC + AC)^ 



with 
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3. pr, -► K*K* 



The amplitude for this reaction is given by 

aw-*- = ^ac* + ac q + ac q + m%t + MD^e 4o 
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4. pr/ -► p/9 



The amplitude for this reaction is given by 
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Jvl 9b 
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(B8) 



5. pr] — > ttuj and pr] — > 7T7T 



The amplitudes .Mpr^™ for the reaction pr? — > nu and .Mp^-^ for the reaction pr] — > 7T7t 
can be obtained from the amplitudes for the reaction 7rr] — ► pc<j and 7rr] — > np via crossing 
symmetry, i.e., interchanging p 4 with either — p 3 or — p 4 . Explicitly, they are given by 



M 
M 



pTj— >7T7V 



Mm^npipi «"> -p 4 )- 



(B9) 



APPENDIX C: ETA ABSORPTION BY OMEGA MESON 

The amplitudes for these reactions can be obtained from corresponding ones for eta 
absorption by p° meson by replacing the mass of the p meson with that of uj meson or by 
the crossing symmetry via interchanging pi with — p 4 , i.e., 



M 



Lorj^KK 



M n 



p°rj-*KK( m p ~ > m uj), 
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ujt)-^KK*(KK*) 
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} (m p -> m w ), 



M 



ujri^K'K* 
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uir)— >irp 



-M pr? _ >7ra; (pi «-> -p 4 ). 



(CI) 



APPENDIX D: ETA ABSORPTION BY K MESON 



The amplitudes for the final states of 7rK, 7rK*, pK, and pK* can be obtained from those 
for the reactions shown in Figs. ^ and El using the crossing symmetry. For the final states 
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of ujK and ujK*, the amplitudes are related to those with final states p°K and p°K* by 
replacing the mass of p meson with that of uj meson, i.e., 



JvX j-c ri >7T K 

" * A I / ' /i A 
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= M wv _ 


>KK*(Pl 
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>kk(Pi <-»■ - 


-Ps), 






-Ps), 


■M-Kti^ujK 


= Mrt,- 
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(Dl) 



APPENDIX E: ETA ABSORPTION BY K* MESON 



Their amplitudes can also be obtained from those for other reactions via the crossing 
symmetry or replacing the mass of p meson with that of uj meson, i.e., 

M K * V ^ttK = Mk V ^ttK*(Pi -P3), 

M K *r,^K* = M^K*K*(P1 <"»• ~Pi), 

M K *r,^pK = M K r,->pK*(Pl ^ ~Pa), 

M K * v ^pK* = M pv ^ K *K*(Pl <~> ~Ps), 

Mk* v ~*ujK = MK*r ) ^p°K('m p — > mj, 

Mk*t,^wk* = M K *r,^p°K*(m p -> m w ). (El) 
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